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/  ABSTRACT 

I 

I 

Because  of  the  recent  revival  of  interest  in  superconductivity,  we 
have  felt  that  a  pedagogical,  yet  concise  review  of  the  fundamental  ideas  of 
the  Bardeen-Cooper-Schrieffer  (BCS)  theory  of  superconductivity  is  timely. 

We  discuss  the  basic  ideas  of  the  BCS  theory  and  list  the  formulae  for  the 
thermodynamic  quantities  of  this  theory  in  the  usual  real-time 
representation.  We  also  represent  these  formulae  in  tjie  alternate 
imaginary- time  representation,  which  makes  ^hands-on1*  numerical  evaluation 
of  the  thermodynamic  quantities  of  interest  readily  accessible.  Finally,  we 
point  out  the  limitations  of  the  BCS  theory  and  describe  some  of  its 
proposed  extensions  for  describing  the  behavior  of  some  of  the  new, 
unconventional  superconductors. 


Accession  For 

NTIS  GRAAI 

o'* 

DTIC  TAB 

D 

Unannounced 

□ 

Justification _ 

By 

'  Distribution/ _ 

f  Availability  Codes 


Dlst 


n 


A 


Avail  end/or 
Special 


The  discovery  of  "high- temperature"  superconductivity  in 


La^  xBa  CuO^  by  J.  Georg  Bednorz  and  K.  Alex  Milller  (1)  in  late  1986  has 
inspired  intense  experimental  and  theoretical  work  worldwide  in  all  areas  of 
superconductivity.  The  study  of  superconductivity  will  therefore,  no  doubt 
become  a  vital  part  of  the  education  of  the  next  generation  of  scientists. 
With  this  perspective  in  mind,  we  shall  attempt  to  provide  a  concise  survey 
of  the  important  ideas  of  the  Bardeen-Cooper-Schrieffer  (BCS)  theory  of 
superconductivity  (2)  as  a  first  step  towards  understanding  more  advanced 
concepts  in  superconductivity  theory  which  have  arisen  with  the  advent  of 
the  "high- temperature"  and  other  unconventional  superconducting  materials. 

In  addition  to  providing  a  cursory  development  of  the  BCS  theory  of 
superconductivity,  we  also  give  a  compilation  of  formulae  which  describe  the 
thermodynamic  properties  of  a  BCS  superconductor.  Our  emphasis  here  is  on 
enabling  the  reader  to  get  to  the  problem-solving  stage  as  directly  as 
possible.  We  refer  the  reader  desiring  a  more  rigorous  treatment  to  the 
various  texts  on  superconductivity  theory  (3).  A  very  useful  resource  of 
references  on  various  aspect  of  superconductivity  can  be  found  in 
Ginzberg  (4) .  We  include  this  report  with  a  discussion  of  the  limitation  of 
the  BCS  theory  and  a  survey  of  the  extensions  of  it  which  have  been  proposed 
to  treat  unconventional  superconducting  materials  such  as  the  heavy- fermion 
and  "high-T^"  superconductors. 

g  02.  THEQEK 

In  this  section  we  briefly  discuss  the  key  points  of  the  BCS  theory 
following  the  treatment  by  Leggett  (5).  The  important  first  step  towards 
formulating  a  theory  of  superconductivity  was  taken  by  Cooper  (6).  He 
considered  the  idealized  case  of  a  pair  of  electrons  in  a  metal,  one  of 
momentum  k  and  spin  t  and  the  other  of  momentum  -k  and  spin  i ,  which  at 
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temperature  T  -  0  are  coupled  by  the  effective  attractive  interaction  g  <  0. 
The  coupled  or  bound  electron  pair  (the  Cooper  pair)  is  of  lower  energy  than 
a  separated  electron  pair,  provided  they  interact  in  the  presence  of  a 
filled  Fermi  sea  composed  of  the  remaining  N-2  electrons.  That  is  the 
momentum  of  the  paired  state,  k,  is  greater  than  the  Fermi  momentum  kp  and 
the  Fermi  sea  blocks  the  electron  pair  from  being  scattered  to  states  below 
kp.  In  analogy  to  the  single-pair  case,  one  might  think  that  an  ensemble  of 
many  such  pairs  might  lower  the  energy  still  further.  Bardeen,  Cooper  and 
Schrieffer  (2)  showed  that  this  is  indeed  the  case;  the  new  state  so  formed 
is  called  the  superconducting  state. 

The  BCS  ground  state  is  usually  referred  to  as  a  spin-singlet,  even- 
parity  s-state.  This  simply  means  that  the  total  spin  S  of  the  pair  with 
individual  spin  components  +h  and  -  h  is  S  -  0,  and  that  the  relative 
orbital  angular  momentum  of  the  pair  is  £  -  0.  Superconducting  states  can 
be  formed  for  other  values  of  i  and,  in  analogy  to  atomic  orbitals,  are 
labeled  as  p-wave,  d-wave  states,  etc.  We  shall  comment  on  these  extensions 
later. 

The  effective  interaction  g  between  the  Cooper  pairs  must  be 
attractive  to  form  the  superconducting  state,  i.e,  the  coupling  energy  must 
be  greater  than  the  Coulomb  repulsion  between  the  pair.  The  mechanism  by 
which  electrons  pair  establishes  the  genre  of  the  superconductor.  BCS 
superconductors  are  characterized  by  an  electron-phonon  mechanism.  An 
electron  passing  through  the  lattice  of  positive  ions  of  a  metal  will  draw 
these  ions  toward  it,  i.e.,  polarize  the  lattice.  Due  to  the  large  mass  of 
the  ions  compared  to  the  electron,  the  lattice  remains  polarized  for  some 
time  after  the  electron  has  pa  sed.  Therefore,  a  second  electron  can  be 
attracted  to  the  net  positive  charge  and  becomes  effectively  coupled  to  the 
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first  electron.  Experimental  measurements  of  quantities  such  as  electrical 
resistivity,  thermal  conductivity  and  the  superconducting  transition 
temperature  can  be  used  to  determine  the  strength  of  the  electron -phonon 
interaction. 

We  shall  now  briefly  review  the  method  of  second  quantization,  which 
will  help  us  in  describing  the  mathematics  of  superconductivity  theory.  A 
non-interacting  N-particle  system  is  characterized  by  wave  functions 
^i(ri,Si) ’  *  **  1*2. •••N,  with  r^  denoting  positions  and  s^  the  spins  of  the 
particles.  The  wave  functions  form  a  complete  orthonormal  basis  set  of 
which  the  states  of  the  system  are  composed.  In  the  second- quantized 
representation  one  does  not  specify  the  wave  functions  of  the  system  of 


particles;  instead,  the  particle  numbers  n^,^...  in  these  states  are 
specified  as  the  system  variables.  The  particle  numbers  are  the  eigenvalues 
of  the  number  operators  n^  given  by  n^  ■  cici’  wliere  c[  *s  t*ie  creation 
operator  for  state  i  and  c^  is  the  corresponding  annihilation  operator.  The 
operator  c^  acting  on  a  state  i  with  n^  particles  in  it,  increases  the 
number  of  particles  to  (n^+1),  whereas  the  operator  c^  acting  on  the  same 
state  decreases  the  number  of  particles  to  (n^-1).  If  we  are  dealing  with 
fermions,  such  as  electrons,  the  creation  and  annihilation  operators  obey 


anticommutation  relations  such  as  c^c^  +  c|c^  -  1,  c^cj+ctc^-  0(i**j),  etc. 


reflecting  the  antisymmetry  of  the  complete  wave  function  of  the  system. 


The  antisymmetry  of  the  total  wave  function  puts  restrictions  on  the 


possible  values  of  the  particle  numbers  n^  (i.e.,  the  Pauli  Exclusion 
Principle) . 


We  shall  treat  the  BCS  theory  at  T  -  0  first  (the  ground-state 


problem) .  Let  us  write  the  wave  function  of  the  system  ♦  as  a  product  of 
pair  states  given  by 


5 


*-n  (uk  +  vk4,clkl)|°>  . 


Here  |0>  is  the  uncorrelated  vacuum  state,  c^c^^  is  the  creation  operator 
for  the  superconducting  pair  state,  is  the  complex  probability  amplitude 
for  the  electron  pair  state  to  be  occupied,  and  is  the  complex 
probability  amplitude  for  the  pair  state  to  be  unoccupied.  The  total 
probability  for  a  state  at  k  must  be  1,  i.e.,  |u^|  +  |v^|  -  1.  Suppose 

t 

initially  that  the  pair  state  (kt,-ki)  is  occupied  and  that  the  state  (kt, 

t 

-ki)  is  unoccupied  where  k'  and  k  differ  by  momentum  q.  The  probability 
amplitude  that  the  system  is  in  the  initial  state  is  v^u^,  •  If  there  is  a 
scattering  process  in  which  the  first  state  becomes  unoccupied  and  the 
second  state  becomes  occupied,  the  probability  amplitude  for  the  transition 
becomes  u^v^.v^u^,  .  Since  eq.  1  implies  that  the  superconducting  state  is  a 
superposition  of  states  containing  different  numbers  of  particles,  we  must 
work  in  the  grand  canonical  ensemble  for  which  the  energy  operator  is 


K  -  H  -  -  (T+V)  -  mN  , 
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c 1  c 

ka  ka’ 


■  N  -  ]  ckacka  and  V  "  "  } 


'k'a  -k'/9  -k0  ka’ 


Here  ^  is  the  chemical  potential,  T  is  the  kinetic  energy,  and  V  is  the 
potential  energy  containing  the  effective  interaction  g,  discussed  earlier. 
It  is  assumed  that  g  is  independent  of  the  frequency  of  vibration  of  the 
ions  of  the  lattice.  The  spin  indices  a  and  /}  can  each  be  t  or  4- ,  and  N  is 
the  total  number  of  electrons  in  the  system  which  is  fixed  by  the  chemical 
potential  /i.  The  characteristic  energy  of  vibration  of  ions  (the  Debye 
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energy  )<Wp)  is  assumed  to  be  much  higher  than  any  excitation  energy  of  the 
electrons  being  considered.  If  the  wave  function  of  the  pair  state  has  the 
form  given  by  eq.  1,  then,  intuitively,  the  average  quantity  <T-/*N>  is 

<T-mN>  -  ^  2  £fc|vk|2  .  (3) 

k 


The  factor  2  arises  because  there  are  two  electrons  in  the  pair.  The 
quantity  is  simply  the  kinetic  energy  of  an  electon  measured  with  respect 
to  the  chemical  potential  (£  —  -  /i) .  To  get  the  average  potential 

energy  <V>,  we  restrict  ourselves  to  taking  averages  of  pairs  of  operators 
(mean- field  approximation)  such  as 


<ct'ct 
kt  -k' 


~  <  W 


>  (ukV 


* 

F  F 
k'rk 


(4) 


The  non- zero  value  of  the  particular  average  given  above  differentiates  the 
superconducting  state  from  the  normal,  metallic  state,  the  averages  of  other 
combinations  of  operator  pairs  being  the  same  as  in  the  normal  state. 

So  finally,  the  average  energy  can  be  written  as 

OO  -  )  2<k|vk|2  •  1st  }  Vk.q  .  (5) 

k  k,q 


It  is  convenient  to  express  u^  and  v^ 
defined  by 


in  terms  of  a  complex  quantity 


VDk  and  \  -  <V<k>/Dk 


(6) 
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with  Ek  ”  +  ^k2,  and  Dk  "  ^l^kl 2  +  (Ek+^k^2,  such  that 


K  I2  +  hJ2  - 1  and  Fk  -  Yk  -  V2Ek  - 


a 

To  find  the  ground-state  energy,  we  minimize  eq.  5  with  respect  to  F^.  By 

noting  that  the  first  term  on  the  RHS  of  eq.  5  can  be  written  as  2  ^  J  A^dF^ 

k 

(apart  from  a  constant)  (5),  the  required  minimization  gives 


\  -  Ul  j  Fk.  "  2jg|  j  Ak,/(2Ek,)  ,  with  k'  -  k-q 


Equation  8  is  the  famous  BCS  self-consistency  relation,  at  T  »  0,  which 
implicitly  defines  the  complex  gap  function  A^.  This  equation,  when  solved 
numerically,  admits  to  a  non- zero  solution  for  the  gap  function  A^.  In  the 
normal  state  the  function  F^  defined  by  eq.  7  is  identically  zero,  implying 
that  A^  -  0.  Thus  A^  serves  as  the  complex  order  parameter  of  the  theory  of 
superconductivity.  For  the  sake  of  simplicity  we  shall  assume  that  A^  is 
independent  of  k  and  drop  its  k- label.  However,  in  the  case  where  the  gap 
function  can  change  with  position,  the  k- label,  as  we  shall  see,  is 
important . 

2 

The  physical  interpretation  of  |a|  ,  introduced  in  eq.  6,  is  that  of 

2  2 

an  energy  gap.  The  relation  E^  -  ,/£*+!  a|  implies  that  there  is  a  gap,  A,  in 
the  energy  spectrum  of  a  superconductor.  This  gap  appears  because  a  single 
pair  cannot  be  broken  without  a  finite  expenditure  of  energy.  Thus  the 
energy  gap  originates  from  the  fact  that  there  is  a  binding  energy  between 
electrons  in  a  Cooper  pair.  The  elementary  excitations  of  the  ground  pair 
(GP)  state  are  of  two  kinds.  The  first  kind  can  be  formed  by  breaking  up 
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the  pair,  i.e. ,  only  one  of  the  two  particles  occupies  the  pair  state. 
Therefore,  there  is  a  single  particle  in  the  kt  or  in  the  -kl  state,  each 
having  kinetic  energy  (relative  to  ,  but  losing  the  pairing  energy 
(i.e.,  F^  -  0) .  This  type  of  excitation  will  be  called  the  broken  pair  (BP) 
state.  The  second  kind  of  excitation,  called  the  excited  pair  (EP)  state, 
can  be  formed  in  a  manner  similar  to  eq.  1,  the  only  difference  being  that 
the  new  probability  coefficients  for  this  state  are  so  chosen  that  this 
state  is  higher  in  energy  and  is  orthogonal  to  the  state  given  by  eq.  1. 

The  GP,  BP  and  EP  states  are  analogous  to  the  bonding,  nonbonding  and 
antibonding  states  of  molecular  hydrogen.  It  can  be  shown  that  the 
excitation  energies  of  a  single  BP  state  (there  are  two  such  states)  and  the 
EP  state  relative  to  the  ground  pair  (GP)  state  are  (for  details  see  Leggett 

(5))  EBP  '  EGP  ”  Ek  and  EEP  ‘  EGP  ”  2Ek  ' 

The  concepts  of  the  three  states,  GP,  BP  and  EP,  help  us  to  generalize 

the  BCS  theory  to  finite  temperatures  (T  >  0)  .  The  quantities  u^,  v^,  A^, 

E^  and  now  become  implicitly  dependent  on  temperature.  The  probabilities 

of  occurrence  of  these  states,  are  given  by 

PGP  “  n-1,  PBp  -  n'1  exp(-^Ek],  P£p  -  n'1  exp[-20Ek]  ,  (9) 

with  n  -  (1  +  2  exp(-/9Ek]  +  exp[-20Z^])  where  0  -  1/T  (using  units  in  which 
Boltzmann's  constant  is  unity).  The  quantity  F^  in  eq.  7  is  generalized  to 

Fk  ”  uifcvk(PGP"PEP^  “  (A/2V  tanh  <*V2)  •  (10) 

Notice  that  the  BP  state  does  not  contribute  to  eq .  10,  since  the  amplitude 
factor  u^v^  i-s  zero  for  this  state. 


The  Real-Time  Representation 

The  evaluation  of  thermodynamic  quantities  within  BCS  theory  can  be 
performed  in  two  ways.  The  first  method,  called  the  real-time 
representation,  follows  along  the  lines  of  the  theory  outlined  in  the 
previous  section  for  T  -  0.  The  second  method,  called  the  imaginary- time 
representation,  will  be  discussed  in  the  next  section.  In  the  real-time 
representation,  the  time  dependence  of  creation  and  annihilation  operators 
involves  exponential  terms  of  the  form  exp(±iKt),  where  K  is  the  energy 
operator  defined  in  eq.  2.  If  the  wave  vectors,  k,  are  closely  spaced, 
forming  a  quasi -continuous  spectrum  near  the  Fermi  surface,  then  we  can 
replace  the  sum  over  k  by  an  integral  over  energy  through  the  correspondence 
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N(0) 


(ID 


where  N(0)  is  the  density  of  states  (i.e.,  number  of  states  per  unit  energy 
interval)  per  single  spin  evaluated  at  the  Fermi  energy  (i.e.,  n-0) ,  and 
is  the  Debye  frequency.  The  self-consistency  relation,  eq.  8,  can  now  be 
generalized  with  the  help  of  eq.  10  to  give 

WD 

A  -  |g|  )  Fk  -  A  J  defc[A/2Ek]  tanh[Ek/2T]  ,  (12) 


where  A  -  N(0)|gj.  At  the  transition  temperature  (T  -  Tc),  A  -  0,  so  that 
eq.  10  reduces  to 
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1  -  A  [  d^k(2|k)'1  tanh(£k/2Tc>  ,  (13) 

'WD 

or  A'1  -  in[rWD(2jrTc)'1]  ,  (14) 

where  inT  is  equal  to  Euler's  constant  (0.577).  Rearranging  eq .  14  yields 
the  transition  temperature  as 

Tc  -  1.135  wD  exp( -1/A)  .  (15) 

At  T  -  0,  A  —  Aq  and  eq.  12  can  be  directly  solved  to  give 
WD 

a'1  -  J  d£  (C2  +  |A0|2)*4  *  in(2u;D/A0)  (v  a>D/AQ  »  1)  .  (16) 

Hence,  Aq  -  2wDexp( -1/A)  (17) 

By  eliminating  A  from  eqs.  15  and  17  we  get  (A^/T^)  -  1.76  which  is  a 
universal  constant  of  the  BCS  theory,  being  independent  of  the  parameters  of 
a  particular  material  such  as  the  density  of  states  and  the  coupling 
constant,  g.  For  temperatures  0  <  T  <  Tc  eq.  12  is  solved  by  eliminating  A 
and  «D  through  the  use  of  eqs.  15  and  17  and  through  the  use  of  the 
transformation  tanh(x/2)  -  1  -  2exp(-x)/[l  +  exp(-x)].  Equation  12  is  thus 
divided  into  two  terms,  the  latter  being  convergent  (i.e.,  integrable  from 
the  limits  0  to  ®) .  We  have  found  that  a  15-point  Gauss -Laguerre 
quadrature  (7)  can  be  readily  applied  to  solve  eq.  12,  and  is  particularly 
accurate  at  low  temperatures  where  w^/T  »  1.  Numerical  solutions  of  eq.  12 


become  quite  inaccurate,  however,  near  T  since  one  has  to  determine  a  small 

c 

quantity,  A,  by  subtracting  two  large  numbers.  The  imaginary- time- 
representation  treatment  avoids  this  problem  and  is  much  simpler  to  use  in 
practice . 

Having  obtained  the  temperature-dependent  gap  function  A(T) ,  as  given 

above,  we  can  now  calculate  the  free  energy  difference  F  -  F  <  0  of  the 

s  n 

superconducting  state  relative  to  the  normal  state.  Noting  that  <V>  - 
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Ffc  ,  we  can  write  (5) 


Fs  -  Fn  -  <V>  +  ^[2i  F*  -  2  J  F*dA]  -  ^<AF*  -  2  J  F*  dA) 

k  k 


-  -  2  N(0)  M2  -  N<°>  |  <^k  xdx(x2  +  ^  )  tanh(7x2+£2/2T)  , 


i  2  r D 

-  -  -  N(0)  | A |  -  2T  N(0)  I  dfk  2n[cosh(Ek/2T)/cosh(£k/2T) ] 


Eq"a*;ion  18  also  gives  the  value  of  the  thermodynamic  critical  field  Hc , 
2 

since  F  -  Fg  -  Hc/(8w).  The  entropy  of  the  superconducting  state  at 
temperature  T  is  given  by 


Ss  “  ^  (PGPinPGP  +  2  PBPinPBP  +  PEPinPEP^ 


u^d 

-  N(0)  |  d^k  ^[1  -  tanh(Ek/2T)]/T  +  2in[l  +  exp(-Ek/T)]|  . 


12 


The  specific  heat  at  constant  volume,  Cs ,  of  the  superconductor  is  easily 
calculated  from  eq.  19  by  the  formula 


C 

s 


dS 

-  T(TT)  -  N(0) 


(Ek/T) 


sech 


(20) 


It  is  useful  to  compare  Cs(T)  above  with  cn(Tc),  the  specific  heat  of  the 

2 

normal  state  at  T  .  Since  C  (T)  -  7T,  where  7  -  2n  N(0)/3  is  the  Sommerfeld 
c  n 

constant,  we  have 


C  (T)  91  9  9  r  n 

-  3(2*  )'  (A/T.)  <T./Tr  2  (-D  (n+l)[K3(xn)  +  3K1(xn)] 

C  n-0 


3(25r2)'1  (i-  ^-)  Y(T)  , 
c 


(21) 


with  Y(T)  -  J”d^k(l/2T)  sech2(Ek/2T)  -  1  -  a 


(22) 


where  x_  -  (n+l)A/T  and  q'L  -  [1  -  (T^~)/(2A2)] 


n 


In  obtaining  eq.  21  we  have  used  the  integral  representation  of  the  modified 

CO 

Bessel  function  K^(x)  and  the  identity  sech2(x/2)  -  4  ^  (-l)n(n+l) 

n-0 

x  exp[ - (n+l)x] .  The  temperature -dependent  quantity  Y(T)  is  called  the 
Yosida  function  and  gives  the  spin  susceptibility  of  the  superconductor. 

The  specific  form  of  Y(T)  in  eq.  22  was  obtained  from  the  discussion  on  the 
electromagnetic  properties  of  the  superconducting  state  in  the  original  BCS 
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paper  (2);  this  term  was  neglected  in  (8).  The  low-  and  high- temperature 
limits  of  these  quantities  are  discussed  in  the  Appendix;  numerical  values 
of  the  thermodynamic  quantities  over  the  entire  temperature  range  are  given 
in  Rickayzen  (9). 

The  Imaginarv-Time  Representation 

In  contrast  to  the  previous  section,  one  can  think  of  temperature  as 
an  "imaginary  time"  r  by  writing  the  Boltzmann  factor  as  exp(-A/T)  - 
exp(ifir)  with  r  -  i/T,  where  r  is  defined  over  a  range  from  0  to  1/T.  The 
variable  conjugate  to  r  is  called  the  Matsubara  frequency,  In  the 

imaginary- time  representation  all  thermal  averages  are  expanded  in  a  Fourier 
series  in  (in  contrast  to  the  Fourier  integrals  of  the  previous  section). 
The  antisymmetry  of  the  fermion  wavefunction  imposes  the  constraint  that 
-  irT(2n+l)  with  n  -  0 ,  ±1 ,  ±2 ,  .  .  .  .  We  can  express  the  formulae  of  the 
previous  section  as  a  series  in  very  elegantly  by  the  method  of  analytic 
continuation  discussed  in  (10). 

Let  f(w)  -  [1  +  exp(w/T)]  ^  be  the  Fermi  distribution  function  for  a 
state  of  energy  w  at  temperature  T,  with  poles  located  at  u  -  tiw^.  Suppose 
we  wish  to  evaluate  the  integral  ji  du  f(w)  F(w)  around  a  contour  which 
encloses  only  the  poles  of  f(w)  but  not  the  poles  of  the  arbitrary  function 
F(w) .  Then  it  can  be  shown  that 

|  tL>  f(w)  F(w)  -  -  2»iT  ^  F(iu>n)  .  (23) 

u> 

n 


As  a  practical  application  of  eq.  23  consider  the  self-consistency  equation, 
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eq.  12.  Since  the  integrand  is  an  even  function  of  we  can  replace 


p  u  p  u 

j  d£^  by  2  J  d£k>  By  changing  the  integration  variable  from  to  one 


D 

obtains 


J0 


/viAi2 


dEk  <VlAr> 


tanh(Ek/2T) 


(24) 


We 


can  write  tanh(Ek/2T)  -  f (E^) [exp(Ek/T)  -  1]  and  use  eq.  23  [note  that 


exp(i«n/T)  -  -1]  to  get 


1  -  2jtAT  ^  [cA|Aj2f4  ,  wn  -  irT(2n+l)  .  (25) 

u>  >0 
n 

At  T  -  T  ,  A  -  0  so  that  we  get  from  eq.  25 

"D 

1  -  A  2  l/(n  +  0.5)  ,  -  NINT[wD(2xTc)'1  -  0.5]  (26) 

n-0 

Here  NINT  stands  for  the  nearest  integral  part  of  the  quantity  inside  the 
square  brackets.  As  n^  -*  «,  the  RHS  of  eq.  26  is  logarithmically  divergent, 
hence  the  need  to  truncate  the  sum  at  the  Debye  number  n^  associated  with 
the  Debye  frequency  u>D  at  Tc .  For  large  n^,  the  series  on  the  RHS  of  eq.  26 
can  be  evaluated  (7)  and  precisely  gives  the  familiar  eq.  14.  In  order  to 
solve  for  A(T)  as  a  function  of  temperature,  we  note  that  eq.  26  can  be 


rewritten  as 
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"d 

A*1  -  in(t)  +  ^  1/ (n+0 . 5)  ,  t  -  T/T,  (27) 

n-0 

Eliminating  A  from  eq.  21  and  27,  we  obtain  the  self-consistency  condition 

"d 

int  +  ^  [ (n+0. 5)  1  -  2wT(u>2  +  | a| 2) -  0  .  (28) 

n-0 

Equation  28  is  clearly  convergent  as  -»  ®,  so  that  the  upper  limit  in  the 
sum  can  be  extended  to  infinity  if  needed.  We  find  that  A(T)  can  be 
calculated  accurate  to  five  decimal  places  with  n^  -  20  near  T£  and  with  n^ 
-  150  at  low  temperatures.  To  solve  eq.  28  numerically,  we  supply  an 
initial  guess  for  A(T)  at  a  given  temperature,  test  whether  it  satisfies  eq. 
25  and  update  it  by  the  method  of  golden  section  (11).  The  procedure  is 
repeated  until  the  desired  accuracy  is  obtained.  The  empirical  relation 

A2(t)  -  A2|l  -  a(t)tr?^t^  exp[b(t)  (t-l)/t]  j-  ,  (29) 

with 

a(t)  -  [1.538  +  1 . 077t2(  t- 1 . 5)  ]  1  ,  r;(t)  -  [2  +  2t2  ( t- 1 . 5)  ]  ‘ 1 

and 

b(t)  -  [0.5669  +  0.2468t3'33(t-1.3)] ‘1  , 


where  t  -  T/T  ,  provides  a  good  initial  guess,  as  well  as  yielding  proper 
values  for  the  gap  and  its  temperature  derivative  at  the  limits  t  -  0  and  t 
-  1. 

In  order  to  find  in  the  imaginary- time  representation,  eq. 

18  is  convenient,  since  the  tanh(E^/2T)  function  can  be  expressed  in  terms 
of  f(E^),  as  before.  Carrying  out  the  analytic  continuation  first  using  eq. 
23  and  then  performing  the  integration  with  respect  to  A,  we  get 


"D 


F  -  F_  -  N(0)  4 -A  in( t)  +  4»r  T 
n  b 


5  "• 


A(2  +  w2)4  -  w_  -  |A|2(2Wn)'1]j- 


n 


n 


n-0 


(30) 


Again,  we  see  that  as  -►  ®  (i.e.  ,  -*■  ®)  ,  the  sum  inside  the  square 

bracket  is  convergent.  The  Matsubara  representation  of  the  entropy  S  and 

the  specific  heat  C  (at  constant  volume)  are  related  to  the  free  energy 

s 

dF  dS 

through  S  -  -  (~p  ,  Cs  -  T(~)  and  can  be  readily  obtained  from  eq.  30  by 

noting  that  both  A  and  depend  on  temperature.  In  practice,  however,  we 

have  found  that  F  -  F  ,  S  and  C  are  most  easily  obtained  by  first 

evaluating  the  free  energy  through  eq.  30  followed  by  single  and  double 

numerical  differentiation  of  F  -  F  to  obtain  the  entropy  and  the  specific 

ns 

heat  (apart  from  the  factors  in  front) . 

The  formulae  given  in  this  section  provide  an  easy  numerical  scheme  to 
calculate  the  gap  and  thermodynamic  quantities  for  a  BCS  superconductor  and 
should  be  useful  for  students  and  teachers  alike.  We  would  also  like  to 
point  out  to  those  who  are  interested  in  more  advanced  topics  that  the 
imaginary- time  method  can  be  extended  to  handle  more  complex  situations, 
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such  as  the  evaluation  of  position-dependent  energy  gaps  (12)  which  occur 
when  a  bulk  superconductor  is  in  contact  with  a  thin  film  of  normal  metal 
(called  a  proximity  junction)  or  in  situations  in  which  an  external  current 
is  passed  through  such  a  junction  (13). 

Limitations  and  Extensions 

The  BCS  theory  is  based  on  a  simple  intuitive  model  for  the  effective 
attractive  interaction  between  electrons,  given  in  eq.  2.  In  that  equation 
it  was  assumed  that  the  interaction  potential  is  attractive  for  wave  vectors 
k  which  lie  within  a  small  spread  of  ±Ak  around  the  Fermi  wave  vector  k^  and 
that  the  potential  is  >  0  otherwise.  In  addition,  the  Debye  energy  is 
assumed  to  be  much  larger  than  the  energy  the  electrons.  This 

approximation  is  called  the  "weak- coupling"  limit,  and  many  superconductors 
fall  into  this  category.  However,  there  are  superconductors  for  which  the 
electron- phonon  coupling  is  strong,  requiring  the  BCS  theory  to  be  modified. 
This  is  done  in  the  "strong-coupling"  theory  of  Eliashberg,  which  utilizes 
the  fact  that  the  electron-phonon  interaction  is  not  instantaneous,  but 
rather  is  delayed  in  time  (retardation).  For  details  of  this  theory  we 
refer  the  reader  to  the  book  by  Mahan  (14).  Other  applications  and 
extensions  of  the  BCS  ideas  can  be  found  in  the  excellent  two-volume  set 
edited  by  R.  D.  Parks  (15).  It  should  be  mentioned  here  that  a  strong¬ 
coupling,  phonon-mediated  model  has  been  proposed  to  explain  the  properties 
of  the  new  "high-Tc"  superconductors  (16).  However,  this  model  deviates 
from  BCS  theory  in  that  the  primary  charge  carriers  are  bosons  (integral 
spin)  rather  than  fermions. 

The  BCS  theory  is  a  microscopic  theory  which  starts  from  the  basic 
idea  of  an  attractive  interaction  between  electrons  in  the  presence  of  a 
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Fermi  sea.  Before  the  microscopic  theory  was  conceived,  macroscopic 

theories  based  on  physical  intuition  were  developed.  The  most  notable  and 

successful  of  these  phenomenological  theories  is  due  to  Ginzburg  and 

Landau  (3).  This  theory  asserts  that  there  exists  a  complex  order  parameter 

2 

0(r)  (not  to  be  confused  with  eq.  1)  whose  square  modulus  |V>(r)|  represents 
the  local  density  of  the  superconducting  electrons.  The  order  parameter 
^(r)  is  zero  in  the  normal  state  (T  >  T  )  and  is  non-zero  in  the 
superconducting  state  (T  <  T  ).  Ginzburg  and  Landau  assumed  that  the  free 
energy  of  a  superconductor  can  be  expanded  in  powers  of  ^(r)  and  its 
gradient  V^(r) ;  symmetry  limits  the  kind  of  terms  in  the  expansion.  It  was 
shown  by  Gorkov  that  the  BCS  theory  of  superconductivity  gives  a  free-energy 
expansion  which  is  equivalent  to  the  Ginzburg -Landau  equation  at 
temperatures  close  to  T  . 

As  mentioned  earlier,  the  BCS  theory  assumes  that  the  bound  electron 

pairs  are  in  a  spin-singlet  (S-0),  even-parity,  s-wave  (i-0)  state. 

Extensions  of  this  theory  have  been  proposed  where  this  limitation  has  been 

lifted.  For  a  triplet  state  (S-l),  the  spins  of  the  electron  pair  have  the 

same  magnetic -moment  projections  (i.e.,  the  spin  wave  function  is 

symmetric) ,  and  antisymmetry  of  the  total  pair  wave  function  requires  that 

parity  of  its  angular  momentum  part  be  odd  (for  example,  i-1)  .  A  triplet  p- 

wave  state  would  exhibit  magnetic  behavior  and  has  a  richer  phase  diagram 

than  an  s-wave  state.  It  has  been  proposed  recently  that  p-wave  pairing 

3 

might  occur  in  "heavy- fermion"  materials  (17).  Liquid  He  exhibits  a 

superfluid  transition  at  extremely  low  temperatures  (at  about  3  mK)  which  is 

3 

associated  with  triplet  p-wave  pairing.  The  theory  of  superfluidity  in  He 
can  be  described  as  an  extension  of  the  ideas  of  Balian  and  Werthamer  (1 4) 
and  is  formally  similar  to  the  BCS  theory. 
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As  another  example  of  the  extension  of  the  ideas  of  the  BCS  theory, 
consider  the  spin-singlet,  d-wave  (2-2)  pair  state,  which  also  has  been 
proposed  for  the  "heavy -fermion"  superconductors.  The  present  authors  have 
developed  a  systematic  symmetry  analysis  to  obtain  a  Ginzburg- Landau  free- 
energy  expansion  for  these  even-parity  superconducting  states  (18,19)  which 
enables  the  calculations  of  thermodynamic  quantities  for  the  heavy- fermion 
and  other  unconventional  superconductors.  A  preliminary  application  of  this 
approach  to  model  the  effect  of  anisotropy  in  the  new  "high- temperature" 
superconductors  has  also  been  pursued  with  some  measure  of  success  (20). 

The  extension  of  the  formal  BCS  theory  to  model  these  materials  has, 
however,  been  met  with  criticism.  Two  major  premises  of  BCS  theory  do  not 
appear  valid  for  these  materials.  First,  the  major  pairing  mechanism  is 
though  to  arise  from  a  magnetic  interaction  rather  than  a  phonon  interaction 
("21-24).  This  effect  could  be  incorporated  into  a  BCS-like  theory  through 
the  introduction  of  an  effective,  magnetic  coupling  constant  gm  >  0  if  it 
were  not  for  the  short  range  of  this  interaction.  It  should  be  recalled 
that  the  BCS-theory  is  a  mean-field  theory,  the  second  major  premise,  and 
consequently  it  ignores  the  local  variations  in  the  magnetic  interaction 
which  appear  to  be  important  in  the  "high-Tc"  superonductors .  However,  in 
spite  of  its  limitations,  BCS  theory  continues  to  be  the  standard  in  the 
field  and  should  be  understood  before  other,  more  elaborate  theories  can  be 
studied. 

VI.  CONCLUSIONS 

In  this  paper  we  have  attempted  to  provide  a  concise  reference  source 
on  the  basic  concepts  and  computational  methods  of  BCS  theory  of 
superconductivity.  It  is  our  hope  that  it  will  provide  the  incentive  for, 
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and  be  helpful  in  the  introduction  of  BCS  theory  into  the  advanced 
undergraduate  or  first-year  graduate  curriculum.  The  material  on 
thermodynamic  properties  are  intended  to  provide  the  reader  with  the 
important  analytical  formulae  describing  the  thermodynamic  properties  of  the 
superconducting  state  as  well  as  supply  a  means  of  access  to  the  more 
advanced  literature.  The  simplicity  of  performing  calculations  in  the 
imaginary- time  formalism  should  be  a  valuable  aid  to  getting  "hands-on” 
experience  in  BCS  theory. 

APPENDIX 

Here  we  present  the  limiting  forms  of  the  thermodynamic  functions  for 
a  superconductor  with  gap  A(T)  at  low  temperature  (T  -*■  0)  and  at  high 
temperature  (T  *  T  ) .  We  also  give  some  additional  mathematical  formulae 
which  serve  as  useful  checks  on  the  numerical  results  obtained  from  real¬ 
time  representation  formulae. 

Case  (1)  T  -  0 

As  T  -»  0,  the  order  parameter  A  »  Aq  and  is  approximately  independent 
of  T.  The  leading  deviation  from  Aq  is  exponential  in  nature  (12), 
indicating  the  existence  of  an  energy  gap  of  the  form 

A2  -  Aq  [1  -  y8*T/AQ  exp(-A0/T)]  .  (Al) 

Similarly  (12), 

Fb  -  F  -  (*2/3)T2N(0)  -  N(0)[^A2  +  exp(-A  /T)  ]  ,  (A2) 

s  n  i  u  o 
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Ss  »  2N(0)(2*AQ/T)h  exp ( -Aq/T) 


(A3) 


and 


Cs(T)/(7Tc)  -  Cs(T)/Cn(Tc) 

-  (6A2)y(x/2)(A0/T)3(A0/Tc)  exp ( -  Aq/T)  ,  (A4) 

with  7  -  2*2N(0)/3. 

The  exponential  term  in  the  specific  heat  is  again  indicative  of  an 
energy  gap  in  the  excitation  spectrum.  Equation  A4  can  be  easily  obtained 
from  eq.  21  by  noting  that  the  temperature  derivatives  are  zero  and  that  we 
can  use  the  asymptotic  form  of  the  modified  Bessel  function  (since  Aq/T  ■* 

00 )  . 


K^(x)  - -  [x/^x)]1*  exp(-x) 

x  -»  <*> 


(A5) 


Case  (2)  T  =  T^ 

Close  to  T  ,  A(T)  is  vanishingly  small  and  serves  as  a  small  expansion 
parameter.  Noting  that  near  T  ,  Y(t)  -*  1  -  2(l-t),  we  get 

A2  *  8»r2l7r(3)  ] _1  T2(l-t)  ,  f(3)  -  1.202  ,  (A6) 

c 

Fs  -  Fn  -  -  4  N(0)  (7f(3)  ]  1  (xTc)2(l-t)2  ,  (A7) 


s  - 

s 


(A8) 


-8N<0)[7f(3)f 1  (*T  )2(l-t)  +  S  , 

c  n 

where  Sn  -  (2jt2/3)N(0)  T, 


and 


c  (T  ) 

s  c 

c  (T  ) 

n  c 


-  1  - 


12. 


7f(3) 


-  1.426 


(A9) 


Note  that  the  RHS  of  eq.  A9  gives  the  jump  in  the  specific  heat  at  T^.  The 
term  -1  on  the  LHS  of  eq.  A9,  which  comes  from  the  nonsuperconductor 
contribution,  is  obtained  by  using  the  asymptotic  form 


Vx) 


id 

2 


A* 

x 


and  the  sum  ^  (-l)n  (n+1)  2  -  w2/12 

n— 0 
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